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We compute the heavy quark potential on configurations generated by the HEM- 



o 

qn^ CGC collaboration with dynamical staggered fermions at 6/^^ = 5.6 and with dynam- 

ical Wilson fermions at 6/^^ = 5.3. The computations are done on 16"^ X 32 lattices, 
corresponding to physical sizes of about 1.6 and 2.3 fm, respectively. Up to the dis- 
tances probed no sign of string breaking is detectable. We also compute the recently 
O proposed scale tq defined by rQ_F(ro) = 1.65. 



The High Energy Monte Carlo Grand Challenge (HEMCGC) collaboration has generated 
ensembles of configurations with staggered sea quarks of mass am = 0.025 and 0.01 at 
6/(/^ = 5.6, and with dynamical Wilson quarks for k = 0.1670 and 0.1675 at Q/g'^ = 5.3, 
on the CM-2 at SCRI, in order to do spectroscopy and to compute some simple matrix 
elements [1, 2, 3]. We have used these ensembles to measure the heavy quark potential in 
the presence of relatively light dynamical quarks. 

Comparing with the heavy quark potential obtained in the quenched approximation, we 
can check for effects from the dynamical quarks. Apart from an interest in its own right, 
such a comparison can help estimate the systematic errors from quenching in recent lattice 
gauge theory computations of heavy quarkonium spectra and the resulting determination of 

Cts- 

The heavy quark potential, or its derivative, the force -F(r), can also be used to obtain 
a scale. As recently suggested by Sommer [4] this scale does not come from the string 
tension, an entity that does not really exist as an asymptotic quantity in the presence of 
dynamical quarks, but rather from the dimensionless quantity rQF[ro) having a fixed given 
value. Sommer chose for this value 1.65, which leads, upon examining the phenomenological 
heavy quark potential, to a scale Tq — 0.5 fm. 

The computation of the potential was fairly standard. We employed the smearing 
method [5] for the space-like parts of Wilson loops W{R^T) to enhance the overlap, Ci, 
with the state with lowest eigenvalue of the transfer matrix in the appropriate channel with 
external charges: 

W{R, T) = ci exp{-ry(i?)} + • • • . (1) 

For the construction of the larger space-like segments of the Wilson loops we also used fuzzing 
or blocking [6]. From the smeared Wilson loops we compute "effective" potentials 

Vt{R) = \og{W{R, T)/W{R, T + 1)) (2) 

If the overlap, Ci, is large the "effective" potentials Vt{R) will become independent of T 
for small T, before the signal is lost in the statistical noise. To illustrate the tremendous 
improvement smearing brings we show in Figure 1 an example of the effective potential from 
planar Wilson loops and from loops with smeared spatial segments as a function of T. As 
can be seen, a plateau appears with smeared loops even at the larger values of i?, but not 
with the normal (unsmeared) Wilson loops. 

We measured smeared Wilson loops, and hence the potential, for all vectors R lying in 
a plane, with the components going up to half the lattice size. The data were taken on the 
HFMCGC configurations with the relevant simulation parameters summarized in Table 1. 
More details about the configurations and the algorithms used to produce them can be found 
in [1, 3]. The potentials are shown in Figure 2. Rotational invariance seems to be restored 
quite well at distances larger than about three lattice spacings. 

As can be seen, in all cases the potential appears to rise linearly at larger distances, up to 
the largest distances probed. It is expected that this linear rise ends (becomes screened) at 
some distance rbreak due to creation of a quark-antiquark pair that will bind to the external 
charges and thus break the string between them. Only the last two points with Wilson 
fermions of k, = 0.1675 hint at such a string breaking, but without statistical significance. 
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Figure 1: The effective potential extracted from planar normal Wilson loops and from, smeared 
Wilson loops. The sample comes from the fi = 5.6, ma. = 0.025 staggered configurations. 
The data are from the top for R = 8 (x) (for which no .signal was found from normal Wilson 
loops), 6 (n), 4 (O) and 2 (o). 
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am or k 


# meas 


AT 


staggered 
staggered 


2 
2 


5.6 
5.6 


0.025 
0.01 


200 
200 


10 
10 


Wilson 
Wilson 


2 
2 


5.3 
5.3 


0.1670 
0.1675 


240 
140 


10 
9 



Table 1: Summary of the configuration samples .studied. AT denotes the number of trajec- 
tories, of unit length, between mea.surements . 
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Figure 2: The potential for the four data sets, obtained from time separations T and T + 1. 
The lines are the best fits, with the solid portion in the range used for the fit. 



fermions 


am or k 


T 


aVo a^a e f 


ro/a roy^ 


staggered 


0.025 
0.01 


4 
5 


0.723( 2) 0.0576( 4) 0.312( 3) 0.38(2) 


4.8(1) 1.15(2) 


0.714( 5) 0.0589( 8) 0.302( 6) 0.35(4) 


4.8(2) 1.16(5) 


4 
5 


0.757( 4) 0.0481( 5) 0.355( 7) 0.36(3) 


5.2(2) 1.14(4) 


0.767( 8) 0.0470(10) 0.373(15) 0.29(5) 


5.2(3) 1.13(7) 


Wilson 


0.1670 
0.1675 


3 

4 


0.785(14) 0.1255(20) 0.366(23) 1.1(2) 


3.2(1) 1.13(4) 


0.800(15) 0.1209(29) 0.378(16) 0.61(5) 


3.2(2) 1.11(7) 


3 

4 


0.812( 5) 0.0976( 9) 0.382( 6) 0.57(4) 


3.7(2) 1.16(6) 


0.805(11) 0.0978(22) 0.374(12) 0.59(6) 


3.6(2) 1.13(6) 



Table 2: Summary of results from, fits to the effective potentials using Eq. (3) on the four 
data ensembles. The last two columns give the scale ro/a determined from rlF[ro) = 1.65 
and the dimensionless quantity ro^/a . 



A rough estimate of the string breaking distance, rbreak, would be rbreakC = 2?7ig, where 
G is the string tension and m^ a hght constituent quark mass, e.g. niq = mp/2. This would 
lead to the estimate rbreak ~ 0.8 fm, which is somewhat shorter than the largest distances 
probed. Alternatively, one can use the results from the light spectroscopy [1, 3] for nip and 
the string tension determined below: this leads to the estimate rbreak/^ ~ 11 and 5.5 for 
staggered and Wilson sea quarks respectively. For Wilson sea quarks we can measure the 
potential for somewhat larger distances than the rough estimate for rbreak and it appears 
that the string breaking sets in, at the earliest, at almost twice the rough estimate. 

Since we do not observe a string breaking we used the same ansatz that is common in 
quenched simulations [7] for fitting the potential 



V{K) = Vo + (jR 



i-K«^(^'-s 



(3) 



where Gl denotes the lattice Coulomb potential. The last term takes account of the lattice 
artefacts present at small distances; it helps in getting good fits that also include rather small 
distances. We used fully correlated fits with the covariance matrix estimated by a bootstrap 
method. In all cases, the best fit values obtained in this way did not differ significantly from 
those of naive, uncorrelated fits. The best fits, selected according to the criterion of having 
maximal "quality", defined as the product of confidence level times the number of degrees 
of freedom divided by the relative error of the string tension, are also shown in Figure 2. 
The fit parameters, in each case for two choices of the distance T from which the "effective" 
potential was taken, are listed in Table 2. 

From the potential we can compute the force between static quarks by taking appropriate 
differences of the potential 

FAf) = V{f)-V{f-e,) (4) 

Recall that we computed the potentials for r in a plane, so we get two non- vanishing com- 



ponents of the force. We are interested in its magnitude F{r) = J F^{f) + F2{r). For the 






0.3 



0.2 



0.1 



0.0 
0.3 



w 0.2 

cd 

0.1 



0.0 



staggered fermions 

^ 5.6, am = 0.025, T = 4 




Wilson fermions 

= 5.3, K = 0.1670, T = 3 




10 



Staggered fermions 

8 = 5.6, am = 0.01, T = 4 




Wilson fermions 

(S = 5.3, K = 0.1675, T = 3 




10 



R/a 



R/a 



Figure 3: The force for the four data sets, obtained from taking differences of the potential 
shown in Figure 2 as described in the text. The line is the derivative of the best fit to the 
potential with f set to zero. 



distance f we choose 



FL"'\f) 



(5) 

where FL^f) is computed as in Eq. (4) but using the lattice Coulomb potential Gl instead 
of V{r). With this choice F{r) is a tree-level improved observable. The force obtained in 
this way is shown as crosses in Figure 3. 

If we know the potential on a ray nJ, we can compute the force along this ray directly 
by 

F{r) = \v{nd) - V{{n - l)d)\ \d\-^ (6) 



where f is again computed from the lattice Coulomb potential. The force computed in this 
way from on-axis potentials is shown as octagons and from the potential along a diagonal 
as squares in Figure 3. The force, obtained from taking differences, is much noisier than the 
potential, and rotational invariance is not obeyed quite as well. 
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a{ro) a{mp) a{mp) 


V^ 


staggered 


5.6 
5.6 
5.6 


0.025 

0.01 

0.0 


0.102( 2) fm 
0.096( 4) fm 
0.092( 7) fm 0.110( 2) fm 0.119( 2) fm 


454( 8) MeV 
450(16) MeV 


Wilson 


5.3 
5.3 
5.3 


0.1670 
0.1675 

He 


0.156( 9) fm 
0.135( 7) fm 
0.116(16) fm 0.109( 2) fm 0.124( 3) fm 


446(16) MeV 
458(24) MeV 



Table 3: The lattice spacing from Vq for the four ensembles, the lattice spacing from. Tq, p 
and nucleon mass in the extrapolation to zero quark mass or Hc, and \/g with the scale set 
by r-Q. 



As suggested by Sommer [4] the force can be used to obtain a scale which determines the 
lattice spacing. He suggested to define Tq from the dimensionless quantity 



rlFiro) = 1.65 



(7) 



and found from phenomenological heavy quark potentials that Tq — 0.5 fm. There are several 
advantages in using Tq to set the scale. It is well defined even in the presence of dynamical 
quarks where an asymptotic string tension does not exist. It does not require any fitting 
but only some interpolation of the force, and hence avoids biases from a particular fitting 
procedure. It comes from a relatively small distance where measurements are usually more 
accurate. The drawback comes from having to take numerical derivatives, which increases 
the noise. 

The values obtained for ro/a and for the dimensionless combination rQ^/a are shown in 
the last two columns of Table 2. The combination rQ^/a does not seem to depend on the 
kind of sea quarks present, or their mass, within about 5% statistical errors. It also seems to 
agree with the results from quenched simulations, compiled in Ref. [8], again within errors. 

We compare, in Table 3, the lattice spacing extracted from the scale Tq = 0.5 fm, with 
the lattice spacing extracted from light hadron masses from the extrapolation to zero quark 
mass [1] or to n^ [3]. We also list the value for the string tension, obtained from using Tq 
to set the scale. For Wilson sea quarks the lattice spacings from Tq and from light hadron 
spectroscopy are in reasonable agreement but the determination from Tq has large errors from 
the extrapolation to Kc- For staggered sea quarks the lattice spacing from Tq is considerably 
smaller than from the light staggered spectroscopy. However, the HEMCGC collaboration 
has also done spectroscopy with Wilson valence quarks on the ensembles with staggered 
dynamical fermions [2] and found lattice spacings, from the rho mass, of a ~ 0.099 fm and 
0.092 fm for quark masses am = 0.025 and 0.01. These values, and their extrapolation to 
am = 0, a ~ 0.088 fm, agree better with the lattice spacings determined from Tq. 

In conclusion, we have measured the heavy quark potential in the presence of two fiavors 
of moderately light sea quarks, both staggered and Wilson. Up to the largest distances 
probed, about 1.5 fm for the Wilson quarks, we did not see any sign of string breaking. 
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Hence this string breaking seems to occur at the earhest at about twice the distance one 
would expect from a naive rough estimate of 0.8 fm. Quahtatively, the potentials look very 
similar to quenched potentials at comparable lattice spacings, as has already be seen in [9] 
for the case of four flavors of staggered sea quarks. Also the dimensionless quantity Tq^/g^ 
where Tq is determined from rQF(ro) = f .65, seems to agree well with quenched results. The 
Coulomb coefficient, e, in our fits, eq. (3), is of comparable magnitude to those obtained in 
quenched calculations. It would be nice to better determine the running coupling constant 
in the presence of dynamical fermions, since there one should see a difference compared to 
quenched calculations [10]. That would require simulations at smaller lattice spacing, to get 
good resolution of the short-distance part of the heavy quark potential. 
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